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We report on a matrix-based diffraction integral that evaluates the focal field of any diffraction- 
limited axisymmetric complex system. This diffraction formula is a generalization of the Debye 
integral applied to apertured focused beams, which may be accommodated to broadband problems. 
Longitudinal chromatic aberration may limit the convenience of the Debye formulation and, ad- 
ditionally, spatial boundaries of validity around the focal point are provided. Fresnel number is 
reformulated in order to guarantee that the focal region is entirely into the region of validity of 
the Debye approximation when the Fresnel number of the focusing geometry largely exceeds unity. 
We have applied the matrix-based Debye integral to several examples. Concretely, we present an 
optical system for beam focusing with strong angular dispersion and free of longitudinal chromatic 
aberration. This simple formalism leaves an open door for analysis and design of focused beams 
with arbitrary angular dispersion. Our results are valid for ultrashort pulsed and polychromatic 
incoherent sources. 

PACS numbers: 03.50.-z, 41.85.Gy, 42.15.Dp, 42.25.Fx, 42.60.Jf 



I. INTRODUCTION 



Interaction of high-intense ultrashort pulsed radia- 
tion with matter enhances multiphoton ionization T, ^ , 
high-harmonic generation 0, 0[, and supercontinuum 
formation]^ 1^. The above experimental demonstrations 
give rise to significant applications such as multiphoton 
excitation fluorescence imaging ^ and high-harmonic 
microscopy ^ [Oj] , where microscopic structures of trans- 
parent samples may be probed in the vicinity of the focus 
of a tightly focused beam. Such broadband radiation has 
intrinsically a serious sensitivity to temporal dispersion 
and chromatic aberrations lo',llH. Consequently, spatio- 
temporal control of broadband light focusing is manda- 
tory. 

Some authors have demonstrated simultaneous spatial 
and temporal focusing of femtosecond pulses by means of 
wavefield division [ij, . This technique is based on the 
spatial separation of the spectral components of pulses 
into a collection of off-axis beams, also called "rainbow 
beam" , thus allowing a parallel processing, and recombin- 
ing these components at the focal point of an achromatic 
objective lens. The proposed arrangements incorporate 
diffr active gratings, frequently combined with refractive 
prisms, with specific dispersive behavior. Similarly, gen- 
eration of multiple spots in the back focal plane of a lens 
results in a notable energy dispersion, and some pro- 
posals with spatial-dispersion compensation properties 
have been presented elsewhere [ij, . Interestingly, the 
diffraction gratings may be appropriately substituted by 
Fresnel lenses for spatial dispersion compensation of in- 
coherent white light 16, 17], profiting from the inherent 
radial symmetry of these optical elements. A complete 
system design may be performed in terms of the ABCD 
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transfer matrix theory ll^. 

The focal region is extensively analyzed in the litera- 
ture under the assumption that the incoming spherical 
wave is monochromatic and the beam is tightly focused 
such that the Fresnel number of the focusing arrangement 
is high . The Fresnel number of the apertured fo- 

cusing lens, N = €? jXf ^ depends on the aperture radius, 
a, the lens focal length, /, and the radiation wavelength, 
A. In a general diffraction problem, the Kirchhoff diffrac- 
tion theory should be used. However, when the Fresnel 
number is much higher than unity, some spatial symme- 
tries in the vicinity of focus are found, and the Debye 
approximation holds 21, 22]. By the way, the concept of 
Fresnel number has recently been interpreted for broad- 
band coherent radiation [2^. However, roughly speaking, 
the Debye diffraction formulation has been ignored in the 
evaluation of the focal field of polychromatic (both ultra- 
short pulsed and temporally incoherent) spherical waves. 
Only a few studies have taken into account the aforemen- 
tioned approximation [2 1 I25 1I23. [27^ . 

In the transfer matrix theory, the Fresnel-Kirchhoff in- 
tegral is expressed in terms of the Collins formula |2^ . 
which proves to be a powerful tool for system analysis and 
design mil^. Obviously, we may include the analysis 
of either monochromatic or polychromatic focal waves 
of complex optical systems. A large number of focusing 
problems encounter the natural field symmetries about 
the focal point and, therefore, they are opened to be ex- 
pressed in terms of the Debye representation. 

The aim of this paper is to investigate the transfer- 
matrix formalism within the Debye approximation, and 
accordingly to present a simplified form of the Collins 
diffraction integral. In particular, we develop the diffrac- 
tion matrix formulation for broadband complex optical 
arrangements. The paper is organized as follows. In 
Sec. II the basic grounds on diffraction of apertured fo- 
cused waves in the Debye approximation are reviewed, 
adding emphasis to polychromatic waves. In Sec. HI the 
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FIG. 1: Schematic illustration of the optical system under 
consideration. For simplicity, we have depicted the input and 
output optical surfaces of the focusing arrangement. 



ray raatrix theory is introduced. In this way, we ob- 
tain the longitudinal and angular dispersion inherent to 
broadband focal waves. Additional constraints among 
the matrix elements are achieved in a telecentric opti- 
cal arrangement where the Debye representation exactly 
evaluates the focal field. In Sec. IV we extend the De- 
bye representation for the evaluation of the focal field in 
diffraction-limited complex optical systems. The region 
of validity of the Debye approximation is derived and the 
role of the Fresnel number to enclose the focal volume is 
pointed out. In Sec. V we analyze some examples. Re- 
fractive and diffractive apertured singlets are studied in 
detail. The dispersive behavior of the numerical aper- 
ture and the Fresnel number is shown. We also discuss 
the dispersion behavior of a version adapted for focusing 
purposes of the optical system of Ref. with nearly 
frequency-independent Fraunhofer patterns. Finally, in 
Sec. VI the main conclusions are outlined. 



II. DEBYE REPRESENTATION OF FOCAL 
WAVES 

Consider a uniform plane wave incident normally onto 
a focusing optical system. In front of the optical system 
we place an aperture, of amplitude transmittance T(ro), 
at a transverse plane hereafter called input plane, which 
governs the amplitude and the extent of the emerging 
beam. Consider a ray of the incident beam travelling 
parallel to the optical axis and passing though the aper- 
ture at a point of the input plane fq = (ro, ^o) (in polar 
coordinates). If the focusing system is aberration- free, 
the emerging ray is driven to the back focal point F lo- 
cated at a distance / from the second principal point R 
(see Fig.^). A reference spherical surface, VK, of centre in 
the focal point F and radius / reproduces the wavefront 
of the emerging focusing wave. When the sine condition 
is satisfied, the trajectories of the incident and emerging 
rays intersect on the reference surface. 

In general, we may describe the direction of prop- 
agation of the emerging geometrical ray by means 
of the unitary three-dimensional (3-D) vector m = 
(sin^cos(p, sin^sin(/:?, cos6>) written in the Cartesian co- 



ordinate system. From the geometry of Fig. Q we deduce 
that the azimuth of the 3-D vector m is given as 
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The sine condition imposes a constraint for the zenith 
giving 
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where g(Q) = sinO is called the ray projection function 
|3lj |. In general, g{0) gives how a ray entering the optical 
system is projected on the wavefront and may have 
different mathematical expressions other than the sine 
function. For example, under the Herschel condition the 
ray projection function is g{0) = 2sin(^/2). However, in 
the paraxial approximation (^ <C 1) the ray projection 
function is unique, g{0) = 0. 

The focal field may be evaluated in terms of the De- 
bye formulation, which considers the interference of plane 
waves with propagation directions given by the vector m. 
The amplitude of the field at a point P of the focal region 
is then written as 
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where Q{6^(p) is the apodization function, r = 
(r cos(/), r sin0, z), and k = co/c is the wavenumber. 
The numerical aperture of the focusing arrangement, 
NA = sin a, is obtained from the equation g{a) = R/ 
where R is the maximum lateral extent of the aperture. 
Whereas |Tp provides the ray density in the transverse 
input plane, the squared modulus of the apodization 
function | Q P gives the ray density over the spherical ref- 
erence surface W. From the consideration of energy bal- 
ance we find a relation of the apodization function 
and the transmittance of the diffracting aperture 
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Then in the focal volume we have a superposition of plane 
waves of amplitude distribution Q and propagating in 
directions given by m. 

When the incident plane wave is polychromatic, we 
may use Eq. ^ to evaluate the amplitude distribution 
in the focal region for the frequencies constituting the 
field spectral range. In general, an optical system is un- 
able to focus the field of different frequencies in the same 
focal point, being distributed along the optical axis (lon- 
gitudinal chromatic aberration). Compensation of this 
spatial dispersion is available, though the focus position 
is not rigorously independent upon frequency. 

Let us write the Debye diffraction integral for a fre- 
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quency uj as follows 
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In the previous equation we have dropped a factor 
exp (ikz), and we use the trigonometric formula 
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Also, we have used Eq. (^) in order to integrate in the 
angular variable ^o- If represents the focal point cor- 
responding to a given reference frequency cjq, the param- 
eter zo{uj) stands for the axial distance from Fq to the 
geometrical focus for the frequency under consideration. 

Moreover, the position of the second principal point H 
is also spatially dispersed, thus affecting to the value of 
the focal length. Even in a case where the focal point is 
strictly independent of the frequency, a dispersive focal 
length f{uj) caused by an axial shift of H would alter 
the angular distribution 6{uj) of plane waves in the fo- 
cal region, as deduced from Eq. This phenomenon 
is denominated as angular dispersion, and is responsi- 
ble of the dispersive character of the numerical aperture. 
Although spatial and angular dispersion are commonly 
neglected in achromatic objectives, we may design some 
optical arrangements with distinctive dispersive behav- 
ior. Finally, the ray projection function g{6) may also 
vary with the wavelength, though it will be neglected in 
the present paper. 

Previously to the use of the Debye diffraction formula 
of Eq. (|5J for broadband beams, the (spatial and angu- 
lar) dispersion attributes of the focusing system should 
be given, represented in the parameters zo{u;) and /(cj), 
together with the ray projection function. According to 
the paraxial approximation, if the numerical aperture of 
the focused wave is sufficiently small, tan and sin may 
be replaced by 0, and the ray projection function is re- 
duced to g{0) = 0. Performing an additional geometri- 
cal transformation = tq/ f in the diffraction integral, 
based on the paraxial approximation of Eq. ^ , the De- 
bye diffraction integral may be rewritten as 



i27rc f 



J f f'- 

J J-TT Jo 



exp 



CO z 



T(ro,(/>o)exp i 



-I — COS(0 - (po) 



UJ_Zo_ 

2cp'\ 
rodrod^o ,(7) 



where we have used Eq. In this case, the 2-D in- 
tegration is performed in the input plane in lieu of the 
reference spherical surface. Importantly, the focal wave- 
field given in the previous integral equation satisfies the 
parabolic wave equation 



III. TRANSFER MATRIX METHOD 

The dispersive behaviour of an aberration-free (zq in- 
dependent of ro for a given frequency) optical system, 
which is characterized by the focal length /(a;) and the 
longitudinal dispersion z^(uj\ is analyzed hereafter for 
focused beams satisfying the paraxial condition (d <C 1). 
Under the paraxial regime, beam propagation in the op- 
tical system may be represented with the ABCD matrix 
method. Thus we follow a procedure to obtain the val- 
ues of / and zq by means of the elements of the ABCD 
matrix of the optical system. The transfer matrix de- 
scribes both the incident and emerging beams in basis of 
geometrical rays. 

Consider a ray incident onto the optical systems that, 
at the input plane, propagates at a height ro with an 
angle Oq with respect to the z-axis. If the ray traverses 
an optical system described with a transfer matrix 
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the emerging ray exits with a height r and an angle 9 
given by 



To 
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When the incident ray propagates parallel to the optical 
axis, ^0 = 0, as we shall assume here on, the exiting ray 
travels with an angle 



= Coro 



(11) 



in the output plane (and any other transverse plane of 
the focal region). This expression should be interpreted 
as the paraxial approximation of Eq. (O . Consequently, 
the focal distance is exclusively given by one element of 
the matrix. 
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Negative values of r and should be interpreted as an 
inversion of 180 degrees with respect to the optical axis, 
accounted by an increment of tt rad in the azimuth angle. 

Since the focal distance is dispersive by nature, /(cj), 
and thus the matrix element Co{oo), a ray emerges from 
the optical system with an angle given by Eq. ([TT|) that 
varies for different frequencies. If we consider a marginal 
ray propagating with a height ro = which denotes 
the maximum radial extent of the diffracting aperture, 
the emerging ray has a maximum angle a, that is simply 
the (paraxial) numerical aperture of the focusing wave (in 
free space). From Eq. ([TT]) we may obtain the dependence 
of the numerical aperture upon frequency. 
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From Eq. (|Tn|) we have that the incident ray propa- 
gating parahel to the optical axis (^o = 0) emerges at 
the output plane with a height r = Aoro- In the future 
we consider the output plane corresponds to the back 
focal plane for a reference frequency cjq, and therefore 
r{uJo) = and Ao{uJo) = 0. Longitudinal chromatic aber- 
ration makes that the position of the focal point varies 
with frequency, what implies that r(a;), and more impor- 
tantly Aq{uj), may differ from zero for a frequency other 
than c^o- Therefore, we should consider the propagation 
from the output plane along a distance zq in order to meet 
the focus. To find the value of zq we consider the matrix 
Mo and a free-space propagation matrix accounted in the 
complete transfer matrix M as 



(14) 



In the back focal plane for a given frequency the value 
of the element A{lj) vanishes. The solution of the linear 
equation A = is 
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Thus, the longitudinal dispersion is determined from the 
elements Aq and Cq of the transfer matrix. In conclusion, 
from Eqs. (|T^ and ([T^ we infer that the spatial (and 
angular) dispersion of the focused beam is exclusively 
characterized by two elements of the 2 x 2-matrix Mq, 
Aq and Co- 
in a general problem explored with the ray matrix 
method, we evaluate the four elements of the matrix Mq. 
However, the determinant gives unity. 



AoDo - BoCo = 1 



(16) 



thus reducing the degrees of freedom. Previously we have 
deduced that we only two elements of the matrix Mq, 
concretely Aq and Co, are necessary to account for the 
dispersive properties of the focused wave. Consequently, 
the Debye representation derived in the previous section 
should impose a further constraint. In this sense, note 
that the scalar diffraction integrals of Eqs. ^ and Q 
give an exact solution of the focused wavefield for tele- 
centric lens systems 32j. In this case, the aperturing 
screen is placed in the front focal plane, what represents 
the input plane. Since —Dq/Cq evaluates the distance 
from the input plane to the front focal plane, a telecen- 
tric optical system satisfies 
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Therefore, such a requirement additionally assures the 
exact validity of the Debye representation. Unfortu- 
nately, wave dispersion frustrates that I^o vanishes for 
a spectral range and, at least, we may impose I^o to be 
sufficiently small, for example, to vanish for the reference 
frequency ujq. 

Alternatively, the Debye approach is commonly exam- 
ined in terms of the Fresnel number, N = Rq/X\Fo\ [21| . 



In our case, Rq denotes the radius of the exit pupil plane, 
which is conjugate of the input plane. Also, Fq represents 
the distance from the back focal plane to the exit pupil 
plane, 

Fo = -§^. (18) 

Finally, the Presnel number of the focusing geometry is 
written as 
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where A = 27tc/u; is the wavelength. In the previous 
equation we have used that the magnification correspond- 
ing to the input (pupil) plane is Dq^. In telecentric opti- 
cal systems, the value of Fq is infinity and, since Do = 0, 
the Fresnel number also tends to infinity. In the next 
section we will show that the Debye representation ac- 
curately reproduces the focused field when the Fresnel 
number reaches values much higher than unity, 1. 

Finally, let us analyze the situation where the longitu- 
dinal dispersion, given in terms of the axial parameter zq, 
is sufficiently small. By using Eq. (|T6j) we may rewrite 
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We have assumed that Ao{uJo) = and I^o is small in the 
neighbourhood of coo- Under these circumstances, zo may 
be expanded into a power series of the term AoDo. In 
the lowest order we find the approximation zo = AoBo, 
which is valid when 

\AoDo\^l, (21) 

that is, for a given frequency uj sufficiently close to the 
reference frequency uJo- Eq. may be compared with 
Eq. (P3|) . what gives the relation 
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which is consistent with the approximation |AoI^o| <^ 1 
together with the constraint of unitary determinant of 
Mo given in Eq. (tTH|) . 



IV. DEBYE APPROXIMATION OF THE 
COLLINS DIFFRACTION INTEGRAL 

A. Matrix-based formalism and region of validity 
of the Debye representation 

In terms of ABCD matrices, the wavefield emerging 
from a diffraction-limited optical system may be deter- 
mined from the Fresnel-Kirchhoff diffraction integral in 
the form of the Collins formula 29], 
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In Eq. (123 7 Uuj{ro) is the amphtude distribution in the 
input plane, where the aperture is placed, which accounts 
for the (occasionally dispersive) complex transmittance 
of the diffracting element and the spectral strength of the 
source. Also, we have omitted a term exp [iuj{L + z)/c], 
where L is the axial distance from the input plane to the 
output plane, the latter being the back focal plane for a 
reference frequency ooq (usually the carrier frequency of 
pulsed beams or the mean frequency of the power spec- 
trum for incoherent white light). Moreover, the spatial 
coordinate z represents an axial distance from the output 
plane, and r± = (r cos (/>, r sin cj)) indicates the transverse 
spatial coordinates. Finally, the matrix Mq is evaluated 
from the input plane to the output plane, which elements 
are, in general, dispersive and therefore depending on 
frequency. The particular selection of the output plane 
guaranties that Ao{uJo) = 0. 

Eqs. {TJ and allow us to evaluate the 3-D ampli- 
tude distribution in the focal region. However, compari- 
son of the paraxial diffraction integral within the Debye 
approximation and the Collins formula manifests some 
differences. First we address our attention to the de- 
pendence upon the axial spatial coordinate of the phase 
terms in the diffraction integral. In the Debye representa- 
tion, the argument of the phase terms varies linearly with 

in opposition to the more complicated dependence in 
the Collins integral. An approach may be performed af- 
ter a first-order series expansion of the argument about 
the point z = 0. Concretely, we may write 

A _ Ao z _ Ao 
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where AqDo — BqCq = 1 of Eq. llT?>jl has been used. The 
physical interpretation of this approximation hes on the 
assumption that the focal field is mostly concentrated in 
the vicinity of focus, satisfying 
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and then we may write B = Bo -\- zDo ~ ^o- More- 
over, in terms of the axial parameter Fq = —Bo/ Do 
given in Eq. (|18|) . which represents the distance from the 
back focal plane to the image of the input plane (exit 
pupil plane), the previous inequality may be simplified 
as 1^1 <C \Fo\. Consequently, Eq. (|26|) involves that the 
focal volume is far from the exit pupil plane. 

We identify a second difference in the presence of a 
quadratic (on r) phase factor in Eq. external to the 
diffraction integral, which is not found in Eq. {TJ. How- 
ever, this term may be neglected under the assumption 
that the evaluation of the focal field is performed in a 
region where {uo / 2c)\D / B\r'^ <C tt, that is. 



The left side of this inequality has the form of a Fres- 
nel number, and next it will be analyzed in detail. In 
resume, a full analogy of the diffraction integrals given 
from the Fresnel-Kirchhoff formulation and the Debye 
representation is found if the analysis of the focal field is 
accomplished in the restricted region given by Eqs. 
and (|2T|l . 

Let us write the Debye diffraction integral in the 
ABCD matrix representation, 

f .U Z o\ f .00 I \ o /^^x 
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In comparison with Eq. Q we find that the focal dis- 
tance is written as / = ^o- From geometrical consid- 
erations we have obtained the same result written in 
Eq. ([T^ . combined with the approximated equation (|22|) . 
From the comparison we also have that zq//^ = Aq/^o, 
that is, zo = AoBo. Again, this result has been previ- 
ously obtained in Eq. under the assumption that 
1^0 1 ^ 1- Alternatively, this inequality may be writ- 
ten as 



\zo/Fo\^l , 
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which is simply the condition given in Eq. (|26|) for the 
specific value of the axial coordinate z = zo. In conse- 
quence, the focus position observed for a given frequency, 
F{uo)^ is to be far from the exit pupil plane, that is, in 
the near-field region of the focused wave. In the case of 
highly-dispersive focusing elements, this condition may 
bring serious restrictions on the spectral bandwidth of 
the incident beam in order to assure the validity of the 
Debye representation. 

We point out that when = 0, the Collins formula 
given in Eq. reduces to the matrix-based diffrac- 
tion integral in the Debye representation of Eq. In 
this case Fo tends to infinity, that is, the exit pupil is 
located infinitely far away from the focal region. As a 
consequence, the validity of the Debye formulation is ex- 
panded into the whole space. In this singular case one 
cannot speak of "Debye approximation" . 

Importantly, we may conclude that the Debye approxi- 
mation imposes bounds on the spatial coordinates where 
the focal wavefield is evaluated, as seen in Eqs. and 
fTf)) . giving an accurate estimation under a limited lon- 
gitudinal dispersion shown in Eq. ((2^ . 



B. Spatial symmetries 

A relevant issue associated with apertured spherical 
waves correctly described in the Debye representation 
is the inherent symmetry about the focal point. From 
Eq. we find that if the wave function in the input 
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plane is real, U^{ro) G 5R, then the focal field is symmet- 
ric with respect to the focal point, 

UUr±.z + AoBo) = -t/:(-r^, -z + Aq^o) , (30) 

where * denotes a complex conjugate. Note that we 
have considered the focal point is found at a distance 
zq = AqBq. a real wave function is encountered, for 
example, when a pulsed unchirped Gaussian beam prop- 
agates through a diffraction-limited optical system with a 
clear circular aperture. In particular. Wolf et al. 21, 2^ 
considered the spatial symmetries of the intensity in the 
focal region of polychromatic spherical waves. They con- 
cluded that when the diffracting screen is purely absorb- 
ing, the focal intensity is centro-symmetric about the fo- 
cus, where the maximum intensity is observed. 

For the sake of simplicity, we follow our discussion un- 
der the hypothesis that the wavefield in the input plane 
may be factorized in the form Uu;{ro) = So{uj)T{rQ). The 
first function represents the spectral strength of the in- 
cident radiation, which depends exclusively on the fre- 
quency. The second is an azimuthally-symmetric func- 
tion of the radial spatial coordinate, which accounts for 
the aperture transmittance. In this diffraction prob- 
lem we assume that the amplitude transmittance of the 
diffracting aperture is frequency independent, what holds 
for purely- absorbing screens. The Collins diffraction for- 
mula in the Debye approximation is finally written as 
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where R is the radius of the aperture. 



C. The role of the Fresnel number 

Previously we have demonstrated that the diffraction 
integral in the Debye representation may be systemat- 
ically employed with a unique spatial limitation (apart 
from the longitudinal dispersion), conceiving a region of 
validity shown in Eqs. (|26|) and ((Tf)) . In principle, we lack 
the guaranty that the region of interest, that is, the focal 
volume, belongs to the region of validity of the Debye 
representation. Here we analyze the conditions an spher- 
ical wavefield, that emerges from an apertured optical 
system of transfer matrix Mq, should satisfy in order to 
enclose the focal volume in the required bounded region. 

For simplicity, let us analyze the case of a diffraction- 
limited optical arrangement with a circular clear aperture 
of radius i.e., T(ro) = 1 Also, in this examina- 

tion we neglect the longitudinal dispersion, and then we 
restrict our analysis to the case where = for a given 
frequency. The field distribution along the optical axis is 
calculated by means of Eq. which particularized to 
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The on-axis intensity \Uuj{z)\^ is maximum at the ori- 
gin, and the first zeros are found at the points of axial 
coordinates 
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and —z\. The highest values of the on-axis intensity are 
found at points comprised between these two zeros, and 
thus we may consider the region of interest is \z\ < z\. 
Validity of the Debye representation in this region implies 
that the inequality of Eq. is satisfied for \z\ < z\^ 
what means that z\ <C |5o/^o|- Equivalently we may 
write 
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that is, 2 <C A^, where N is the Fresnel number given 
in Eq. (|T^ . In accordance we conclude that a Fresnel 
number much higher than unity guaranties that the focal 
field is accurately described with the Debye representa- 
tion along the optical axis and neighbouring points. 

On regards the off-axis points of the focal volume, let 
us investigate the wavefield in the transverse focal plane, 
z = 0. From Eq. ^^\\ we finally have 
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where Ji is a Bessel function of the first kind and order 1. 
The field distribution given in Eq. is simply the so- 
called Airy disk. Again, the maximum of intensity occurs 
at the origin, and the first zero of intensity is found at 
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Note also that ri gives the limit of resolution of a 
diffraction-limited imaging system in agreement with the 
Rayleigh criterion. Again, the Debye representation is 
valid if Eq. ((TT)) is fully satisfied for points r < ri. Thus, 
we should impose that <C X\Bo/Dq\, what may be 
written as 0.37 <^ N. As a result, high values of the 
Fresnel number allow the focal field to be confined in the 
region of validity of the Debye representation, satisfying 
simultaneously the inequalities of Eqs. and 



V. EXAMPLES 

A. Apertured thin lens 

Let us first consider a positive thin lens of focal dis- 
tance f{uj). The term /o is assigned to the value of the fo- 
cal distance for the reference frequency, i.e., /o = f{^o)- 
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Output 
plane 




FIG. 2: Diagram illustrating an apertured thin lens. The 
diffracting aperture is placed at a distance d in front of the 
lens. The lens focal length is /o for the reference frequency 



TABLE I: Sellmeier coefficients for a common borosilicate 
crown (BK7) glass 34]. 



Coefficient 

Bl 

B2 

B3 

CI 

C2 

C3 



Value 
1.03961212 
2.31792344 10"^ 
1.01046945 
6.00069867 10"^ 
2.00179144 10"^ 



1.03560653 10^ /xm^ 



A diffracting screen is placed in front of the lens at a 
distance as seen in Fig. [21 A negative value of d may 
consider an aperture located at the back of the lens. As 
previously established, the input plane is set in the place 
of the diffracting aperture, whereas the output plane cor- 
responds to the back focal plane for uoq. In this case, the 
transfer matrix Mq has the following elements: 



= 1 



/o 



) = fo^d- 



1 



dfo 



(37a) 
(37b) 
(37c) 



and, obviously, Co{u;) = —f~^{uo). Supposing that the 
focal distance is independent upon frequency, what oc- 
curs for achromatic objectives, the term Aq takes val- 
ues identically zero. However, Aq may reach significant 
values when the focusing lens is highly dispersive, thus 
producing a notable longitudinal dispersion. Also, when 
the diffracting aperture is at the front focal plane for a 
given frequency cj, being d = /(cc;), the focusing system 
is telecentric and, correspondingly, Dq = 0. 

Two sort of thin lenses are analyzed in this sec- 
tion: kinoform-type diffractive lenses and refractive thin 
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FIG. 3: Wavelength dependence of the product Ao^o in a 
central band of the visible spectrum. The elements of the 
transfer matrix are evaluated for a Fresnel lens (solid line) 
and a BK7 glass lens (dashed line) with the same focal length 
/o = 10 cm for a wavelength Ao = 589.3 nm. The aperture 
is placed at: (a) d = /o, and (b) d = 0. The values for the 
Fresnel lens have been multiplied by a factor of 10~^. 



lenses. Diffractive singlets are zone plates which may 
achieve a high-efficiency performance with phase-only 
multilevels 35]. Recently, some authors have addressed 
special attention to investigate the spatio-temporal re- 
sponse of broadband ultra-short pulses focused with zone 
plates 23, 36, 37]. Here we are giving a simplified formal- 
ism to evaluate the focal field in the frame of the Debye 
representation. The focal length of a Fresnel lens may be 
modelized as |33 



/o 



(38) 



Secondly, quartz and glass thin lenses are considered. 
Piano convex lenses are the most basic optical elements. 
They have positive focal lengths and close to the op- 
timum shape for use as focusing lenses for collimated 
beams. If n(uj) denotes the refractive index of the trans- 
parent material the lens is made of, the focal length of 
the refractive singlet may be expressed as 



/o n{uj) - 1 



no 



1 



(39) 



being no = n(cc;o). Bi-convex and other spherical lenses 



8 



may be treated similarly. Instead of using tabulated data 
of the refractive index n{u;) of a given medium, we may 
employ the Sellmeier equation, which is an empirical re- 
lationship between refractive index n and wavelength A. 
The usual form of the equation for glasses is: 



n2(A) = 1 + 



A2-Ci A2-C2 A2-C3 



(40) 



where Bi and Ci, for z = 1, 2, 3, are experimentally deter- 
mined Sellmeier coefficients. Note that this A = 27rc/iJ 
is the vacuum wavelength. As an example, here we con- 
sider a common borosilicate crown glass known as BK7, 
and the Sellmeier coefficients are shown in Tabled 

In Fig. [31 we investigate the longitudinal chromatic 
aberration in relation to proximity to the exit pupil plane. 
This point is important since the Debye approximation 
is restricted to a limited axial dispersion expressed as 
1^0 1 <C 1, as given in Eq. ((2^ . In the figure we 
use the reference frequency coq associated to the aver- 
age wavelength Aq = 589.3 nm of the Sodium doublet. 
In general, the highly-dispersive character of the Fresnel 
lens severely restricts the spectral range we may use the 
Debye representation of the focal field. A relevant prop- 
erty is that when d = /q, i.e. the system is telecentric for 
the radiation frequency cjq, the term AqDq is stationary. 
In other words, AqDq = for the reference frequency loq 
(since Ao{uJo) = 0) and 



duo 







(41) 



for uoq. Consequently, the values of AqDq remain signif- 
icantly low in the spectral neighbourhood around cjq, as 
seen in Fig. |3fa). However, when the aperture is at the 
lens plane ((i = 0), uoq is not a stationary point of AqDq^ 
whose values may increase up to one order of magnitude, 
as shown in Fig. |2fb). The analysis of a BK7 glass thin 
lens has also been included in the figure. Comparatively, 
the refractive lens gives values of IAqI^oI considerably 
lower than those provided by a Fresnel lens. In this case, 
application of the Debye approximation is less stringent 
in terms of bandwidth, even considering the case of aper- 
tures located at the lens plane. 

Validity of the Debye approximation is also restricted 
to Fresnel numbers much higher than 1, what is under 
study in the following example. Consider a circular clear 
aperture of diameter 2R = 1 cm and a thin lens of focal 
length /o = 10 cm for the reference wavelength Aq = 
589.3 nm. If the diffracting aperture is placed at the lens 
plane, the Fresnel number for the reference frequency 



(42) 



In the present numerical example A^o = 424. For other 
frequencies and positions of the aperture, the Fresnel 
number N{u;) varies according to Eq. (|19|) . Let us exam- 
ine a reduced Fresnel number, normalized to the Fresnel 
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FIG. 4: Spectral dependence of the reduced Fresnel number 
N/No when the diffracting screen is placed at (a) d = /o, and 
(b) d = 0. The elements of the transfer matrix are evaluated 
for a Fresnel lens (solid line) and a BK7 glass lens (dashed 
line). The values for the Fresnel lens have been multiplied by 
10~^ in subfigure (a). 



number A'o, written as 



N 



fo 



1 



uo \Bo{u;)\ \Do{uj)\ 



(43) 



When the aperture is placed at the lens plane, d = 0^ the 
terms Bq = fo and I^o = 1 are spectrally invariant (non- 
dispersive), and the reduced Fresnel number is finally ex- 
pressed as N/Nq = uo/uoq. This result is valid for any thin 
(Fresnel and dispersive glass) lens, as seen in Fig. ^b). 
In the case the focusing system is telecentric (for cjq), 
where d = /o, the value of Do{uJo) vanishes and, con- 
sequently, the Fresnel number diverges for the reference 
frequency. As shown in Fig.^a), the Fresnel number for 
neighboring frequencies reaches extremely high values. 
Comparatively, the Fresnel number associated with the 
Fresnel lens are several orders of magnitude higher that 
those of a BK7 glass lens. In general, we may conclude 
that the requirement A^ 1 is fully satisfied within the 
visible spectral band. 

The angular dispersion of the focal waves produced by 
a thin lens may be examined with the matrix element 
Co{uj) and Eq. (|T^ . For instance, a diffractive lens has 
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FIG. 5: Wavelength dependence of the relative numerical 
aperture a/ao for: (a) a Fresnel lens and (b) a BK7 glass 
lens. The curves are evaluated by means of Eq. (CSI (sohd 
line), and Eq. for (i = (dashed line) and d = fo (dotted 
line) . 



high-dispersive numerical aperture given by 



ao 



CO 



(44) 



irrespective of aperture position. Alternatively we may 
use that BqCq ^ —1 in the Debye approximation, and 
thus the spectral dependence of the numerical aperture 
may be expressed as 



(45) 



In Fig. O we plot the relative numerical aperture for a 
diffractive lens (subfigure a) and a BK7 glass lens (sub- 
figure b). The exact expression given in Eq. shows 
that the relative numerical aperture, and thus the an- 
gular dispersion, is independent of the variable i.e., 
the aperture place. However, Eq. ((l5j) reveals a small 
dependence on (i, mainly in the case of the zone plate. 
This point is important since we deduce from Eq. 
that, in the case (i = 0, the numerical aperture is nondis- 
persive, i.e., a{uj) = a^. Strictly speaking, this result 
is false though the angular dispersion is relatively small 
and may be ignored. On the contrary, the approximation 
given in Eq. is extremely accurate when d = /q. We 



may conclude that deviations of Eqs. (|T^ and in- 
trinsically elucidate the degree of accuracy of the Debye 
approximation. 

In Fig. ini we show some contour plots of the three- 
dimensional intensity distribution of an apertured Fres- 
nel thin lens and a BK7 glass lens of focal distance 
/o = 10 cm for a wavelength Aq = 589.3 nm. The 
field intensity is numerically evaluated from Eq. for 
T(ro) = 1 and different frequencies. Also, we have se- 
lected a value of such that the intensity at the focus 
is unity, and the gray levels are plotted in a logarithmic 
scale of intensities to increase the contrast. Again, the 
diffracting screen is a clear circular aperture of diame- 
ter 2R = 1 cm and placed at a distance d = fo. The 
predicted symmetries of the intensity around the geo- 
metrical focus are evident in the plots. In general, we 
observe that the glass lens presents a higher resistance 
to dispersion-induced focal shifts. Additionally, for in- 
creasing frequencies, the focal shift is driven toward the 
lens for the refracting singlet, in opposition to the case of 
the Fresnel lens. Finally, an additional spatial dispersion 
effect is distinct in the subfigures of the glass lens. In or- 
der to have comparable longitudinal dispersion with the 
Fresnel lens, the selected frequencies considerably differ, 
with 20 per cent of increment {Aco/loq — 0.2). In this 
case, the size of the focal spot is notably higher for lower 
frequencies, what is clearly observed by computing the 
number of minima (black shades) along the optical axis 
and in the transverse focal plane. 



B. White-light focusing system with compensated 
spatial dispersion 

A number of designs of achromatic doublet lenses 
[3^ and hybrid refractive-diffractive achromats 40] have 
been proposed to get rid of longitudinal chromatic aber- 
ration inherent in singlet lenses. Ray matrices of this sort 
of focusing setups have (at least approximately) a vanish- 
ing term A{^(uj) = 0, a performance that can be achieved 
over a broad spectral band. Commonly, this property is 
accompanied with spectral invariance of the focal length, 
meaning that both and Co are also invariant in the 
considered spectral range. When used on-axis, an achro- 
matic lens focuses a parallel input beam to a single point. 
However, the off-axis performance is significantly worse 
than the on-axis performance since it is limited by the 
effects of diffraction. For instance, the limit of resolution 
T\ given in Eq. and frequency are inversely propor- 
tional. 

In Ref. 18] Lancis et al. proposed an optical sys- 
tem composed of a nondispersive objective coupled with 
a pair of kinoform-type diffractive lenses to overcome 
longitudinal chromatic aberration and, additionally, the 
diffraction-induced spatial dispersion of the Fraunhofer 
pattern observed in the transverse focal plane. A first 
diffractive lens (DLi) of focal length Zi (for uoq) is placed 
in front of the achromat (of focal length fo for all fre- 
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FIG. 6: Contour plot of the field strength in log scale (natural logarithm of \Uuj{r)\'^) in the focal region of (a-c) a Fresnel lens 
and (d-f) a BK7 glass lens, both of focal length fo = 10 cm and with an aperture placed at d = /o. We consider a frequency 
(a, d) (1 — A)a;o, (b, e) ouq, and (c, f) (1 + A)a;o. In this case A = 0.01 for the Fresnel lens, and A = 0.2 for the glass lens. The 
contour interval (difference between successive contour lines) is 0.75. 



quencies), at a distance s = a -\- and the second 
diffractive lens (DL2) of focal length Z2 is located be- 
hind, at a distance s' = fo -\- o! as depicted in Fig. d 
Under the conditions + = fo^ (lens formula) 
and Z1Z2 = -{s' /sf, the emerging beam is focused at 
a distance /o from the achromatic objective, specifying 
the output plane. In this case, Aq{uo) = in the spectral 
domain. Additionally, inserting the diffracting aperture 
at a distance 



a ■■ 



2 



(46) 



from the first diffracting lens, the rest of matrix elements 
are 

B,{oo) = /o (1 - , (47a) 

^ d ^ Zi^2{d-fo) uo 
fo 4/0 CJ 



(47b) 



and, obviously, Cq = -B^^ (see Eq. ^^). In Fig.Qa' 
is positive and, therefore, the focus is virtual. However, 
inclusion of imaging lenses in the rear may generate a 
real focal plane. 



Eq. (|47a|) indicates that the transfer matrix of the pro- 
posed system have a highly-dispersive element Bo{u;). 
Fig. IHl depicts the spectral dependence of the term Bq^, 
that is, the angular dispersion of the focused waves as 
deduced from Eq. (|45|) . We point out that Eq. (|45|) eval- 
uates the paraxial numerical aperture exactly since Aq 
vanishes. For frequencies sufficiently close to coq we find 
that 



CO 



(48) 



as occurring with a single zone plate. In other words, the 
numerical aperture depends approximately linearly upon 
wavelength. However, whereas a zone plate achieves this 
characteristic angular dispersion at the cost of longitu- 
dinal chromatic aberration, the proposed system has a 
frequency-independent focal plane. As demonstrated in 
Ref. 18], this fact is a direct consequence of imposing 
the reference frequency ujq to be a stationary point for 
the term co/Bq^ that is. 



d_ 

duo 







(49) 



for uoq. This allows a first-order achromatization of 
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FIG. 7: Optical system with compensated spatial dispersion. 
The diffracting aperture is placed at a distance d in front of 
the achromatic objective of focal length /o. The appropriate 
insertion of a pair of kinoform-type diffractive lenses (DLi 
and DL2) is able to strongly alter the angular dispersion of 
the focused beam without introducing longitudinal chromatic 
aberration. 
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FIG. 8: Wavelength dependence of the relative numerical 
aperture for the system in Fig. |3 The thin solid line cor- 
responds to a system with perfectly-compensated spatial dis- 
persion. In the latter case, the numerical aperture has a linear 
dependence with the wavelength, a/ao = X/Xq. In the plot 
we assume Ao — 589.3 nm. 



the Fraunhofer diffraction pattern. Finally, when the 
diffracting screen is a clear aperture, the spectral com- 
ponents of the field (Airy disk) has a nearly invariant 
size and, according to Eq. ((36j) . a frequency-independent 
limit of resolution ri. 

The angular dispersion of this singular optical arrange- 
ment is independent of the parameter (i, that is, the dis- 
tance from the aperture to the achromatic lens may be 
arbitrarily changed. However, the Fresnel number may 
vary substantially for different values of d. The Fres- 
nel number given in Eqs. and is plotted versus 
wavelength for d = fo and d = fo-\- Zi/2 in Fig.lHl In the 
first case (plotted in solid line), the optical system has a 
finite Fresnel number for the reference frequency ujq and, 
therefore, locating the aperture at the front focal plane of 
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FIG. 9: Spectral dependence of the Fresnel number N{oo), 
normalized to the factor A^o(/o/^i), when the diffracting 
screen is placed at d = fo (soUd line), and d = /o + Zi/2 
(dashed line). For the sake of clarity, the values correspond- 
ing to the solid line have been multiplied by 10^. 



the objective does not guarantee a telecentric lens design. 
Contrarily, the second case (plotted in dashed line) ac- 
counts for the condition of telecentricity for cjq, for which 
the Fresnel number diverges. This favours that the spa- 
tial constraints of the focal volume, associated with the 
Debye approximation, are fulfilled. 

Finally, Fig. [1171 depicts the isophotes in the neighbour- 
hood of the focal point for different frequencies. In all 
cases, the Airy disk formed in the focal plane has a com- 
parable magnitude. However, an increase in frequency 
induces a broadening of the focal spot along the optical 
axis. We point out that this behavior substantially differs 
from a conventional achromat lens, where transverse and 
axial widths of the central lobe have a linear dependence 
with the inverse of the frequency. 



VI. CONCLUSIONS 

We have derived an ABCD matrix-based generaliza- 
tion of the Debye representation of focal fields. This 
diffraction integral has been achieved using two alterna- 
tive procedures. In the first approach we substitute the 
focal length and the longitudinal dispersion, which are 
explicit parameters in the classical Debye integral, by 
equivalent expressions given in terms of the ray-matrix 
elements. The second method consists in applying sim- 
ple approximations into the Collins diffraction formula. 
The latter one has the advantage of providing some an- 
alytical constraints which determine the validity of the 
proposed integral. In this sense, the evaluation of the 
focal field is restricted to a bounded region around the 
focus (for a given frequency). As in the standard anal- 
ysis, the inclusion of the focal region into the region of 
validity of the Debye approximation is guaranteed when 
the Fresnel number of the focusing geometry greatly ex- 
ceeds unity. Under the above assumptions, diffraction 
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FIG. 10: Contour plot of the field strength in the focal region of the optical system in Fig. |3 The contour intervals vary 0.75 
in a log scale to augment the contrast. The achromatic objective has a focal length / = 10 cm, and the aperture is placed at a 
distance d = 10 cm. Each plot corresponds to the frequency (a) (1 — A)a;o, (b) cjq, and (c) (1 + A)a;o, with A = 0.2. 



of purely-absorbing screens generates focal fields whose 
spectral components show symmetrical amplitude distri- 
butions around the focus. 

In a second stage, the matrix-based Debye integral has 
been applied to investigate the spectral response of an 
apertured thin lens, either diffracting or refractive. The 
position of the pupil aperture has a strong relevance in 
the Fresnel number, being the telecentric geometry the 
most appropriate design for the Debye representation. 
Due to the longitudinal chromatic aberration, an angular 
dispersion effect is observable and may be important in 
the case of Fresnel lenses. 

Going one step further, we have also examined the op- 
tical setup proposed by Lancis et al. in the frame of 



the Debye representation. This system shows a highly- 
dispersive numerical aperture, approaching the angu- 
lar dispersion of an apertured Fresnel lens, but being 
free of longitudinal chromatic aberration. The numeri- 
cal aperture of the focusing system decreases with fre- 
quency, approximately following an inverse dependence. 
The dispersive nature of diffraction is then compensated 
in the transverse focal plane, showing nearly frequency- 
independent Fraunhofer patterns. As a collateral effect, 
a positive magnification of the focused field along the 
optical axis may be recognized. 

In short, the present formalism leaves door open to 
the analysis and design of focused beams with variable 
angular dispersion. 
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